On the mechanism of branching in negative ionization fronts 
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Abstract 

We explain a mechanism for branching of a planar negative front. Branching occurs as the 
result of a balance between the destabilizing effect of impact ionization and the stabilizing effect 
of electron diffusion on ionization fronts. The dispersion relation for transversal perturbation is 
obtained analytically and reads: s = \k\ / {2(1 + \k\)] — D\k\ 2 , where D, which is assumed to be small, 
is the ratio between the electron diffusion coefficient and the intensity of the externally imposed 
electric field. We estimate the spacing A between streamers in a planar discharge and deduce a 
scaling law A ~ D 1 / 3 . 

PACS numbers: 52.80.Hc, 05.45.-a, 47.54.+r, 51.50.+V 
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One of the greatest unsolved problems in the physics of electric discharges is the clari- 
fication of the mechanism of branching. When a strong electric field is applied to a non- 
conducting medium such as a gas, narrow channels of ionized matter called streamers may 
be formed. The phenomena has been observed in a wide range of scales, from small-gap 
discharges of few centimetres up to scales of kilometres such as in sprites discharges from a 
thunder cloud. The pattern of this branching resembles the ones observed in the propaga- 
tion of cracks, dendritic growth and viscous fingering. Those phenomena are known to be 
governed by deterministic equations rather than by stochastic events. In this paper, for the 
first time, it is derived a quantitative prediction of the branching characteristic length based 
on a deterministic model. The results give an explicit dependence of branching with the 
electric field and the pressure of the gas as it has been observed qualitatively in experiments. 

We use a fluid approximation to describe electric breakdown in non-attaching gases such 
as nitrogen. In these gases, there are indications showing that the most important source 



of electrons and positive ions from neutral molecules is impact ionization Q]. This is the 
only process that will be taken into account in this work. This leads to a minimal streamer 
model from which the electron N e and positive ion N p densities can be obtained. The 
minimal model reads 

° X V ■ (N e S + D e ViV e ) + NelSle'VW, (1) 



Ot 



^ = iV e |£|e-VI*l (2) 
N p - N e = V • E, (3) 

in which S is the local electric field in the gas and D e is the electronic diffusion coefficient. 
The first equation means that the electron density varies in a point as a result of (i) the 
electric current contribution V(iV e £), (ii) the electronic diffusion term D e V 2 N e , and (iii) 
the impact ionization source term A^ e |£|e _1 /' £ '. The second equation means that the positive 
ion density in a point of the gas only varies in time as a result of impact ionization, since the 
mobility of ions is much smaller than that of the electrons. The third equation is Poisson's 
law for the electric field, which is supposed to be irrotational since the magnetic effects are 
neglected. 

All the quantities in Eqs. ([I}-© are dimensionless. The scales of these quantities for ni- 
trogen [2I depend on the gas pressure p. The characteristic field is £ = 200 kV/ cm (p/1 bar). 
The characteristic length is Rq = 2.3 /im (p/1 bar) -1 . The characteristic velocity is Uq = 



76 x 10 6 cm/s. The characteristic time is tq = 3ps(p/lbar) 1 . The characteristic par- 
ticle density is Nq = 4.7 x 10 14 cm~ 3 bar) 2 . The characteristic diffusion coefficient is 
D = 1.8 x 10 4 cm 2 /s(p/lbar)- 1 . 

In the minimal model, a spontaneous branching of negative streamers has been observed 
numerically |3j, as it occurs in experimental situations |4|. In order to understand this 
branching, the dispersion relation for transversal Fourier-modes of planar negative shock 
fronts (without diffusion) has been derived j^. For perturbations of small wave number k, 
the planar shock front becomes unstable with a linear growth rate proportional to k, but all 
the modes with large enough wave number k seem to grow at the same rate. 

In this paper, we address the problem of the selection of a particular wave number in the 
perturbation of negative planar fronts. We will obtain a new dispersion relation depending 
explicitly on the electric field and the electronic diffusion coefficient (which depends also 
on the gas pressure). Our analysis will show that the electron density N e may develop 
steep fronts of thickness 0{\/TT e ) ) satisfying an equation analogous to Fisher equation A 
surprising fact established during the last 30 years is that the combination of sharp interfaces 
with small diffusive effects may result in asymptotic limits (for D e <C 1) in which the motion 
of the interface is described by equations involving solely geometrical properties such as its 



mean curvature ML. These analysis concerned a model today known as Allen-Cahn equation. 
Subsequent work |8j showed that the points of the interface separating two different phases 
move along the normal direction with a velocity proportional to its mean curvature. This 
kind of dynamics is termed " mean curvature flow" . Many mathematicians have contributed 
to provide a rigo rous proof of the convergence of Allen-Cahn model to motion by mean 
curvature 0, Remarkably, some of these limiting models may be such that the 

solutions develop branch-like patterns. In this work, we exploit some of these ideas to study 
the motion of ionization fronts. We will show that a planar front separating a (partly) 
ionized region from a region without charge is such that small geometrical perturbations 
in the charge distribution lead to a motion of the front affected by two opposed effects: 
electrostatic repulsion of electrons and electron diffusion. The first effect tends to destabilize 
the front while the second acts effectively as a mean curvature contribution to the velocity of 
the front, thus stabilizing it. The net result is the appearance of fingers with a characteristic 
thickness determined by the balance of these two opposing actions. 

In order to study the evolution and branching of ionization fronts, we consider the fol- 



3 



lowing experimental situation. The space between two large planar plates, situated at x = 
and x = d respectively (x is the vertical axis and we suppose that d > 1), is filled with 
a non-attaching gas like nitrogen. A stationary electric potential difference is applied to 
these plates, so that an electric field is produced in the gas. The initial electric field is 
directed from the anode to the cathode, along the negative x axis, and is uniform in the 
space between the plates. To initiate the avalanche, an initial seed of ionization is set near 
the cathode. We first study the evolution of planar negative ionization fronts towards the 
anode. 

We will concentrate in the study of the dynamics under the effect of strong external 

electric fields. We denote the modulus of the dimensionless electric field at large distance 

from the cathode as Soo and we will assume that ^> 1. Under these circumstances, 

it is natural to rescale the dimensionless quantities in the minimal model as £ = E, 

N e = £oon e , N p = SooUp, and r = t/S^. For ^> 1, this system can be approximated by 

on 

-^-V-KE + DVHe) = n e |E|, (4) 

V ■ E = n p - n e , (6) 

where D = De/Soo is, in general, a small parameter. Our approximation will be valid in 
all regions where |E| = 0(1). These are the regions of interest since the electric field is not 
expected to vary much in the neighbourhood of the ionization front and we will show that 
it is in this region where the mechanisms leading to branching take place. 

In the evolution of the ionization wave along the x axis, the rescaled electric field can 
be written as E = Eu x where E < 0, so that |E| = \E\ = —E. It is very simple to 
compute travelling wave solutions when D = 0. It can be shown 12j that these solutions 
exist for any c > 1. We will be interested in the limit c — > 1 since it is well known that a 
compactly supported initial data (representing a seed of ionization located in some region) 
develops fronts travelling with this velocity. In the case c = 1 the solution can be obtained 
straightforwardly, giving (with £ = x — ct and c = 1) 

. -e«, for £ < f 1, for £ < { 1 - at, for £ < 

0, for £ > 0, for £ > 



-1, for £ > 
In the case < D <C 1, it is known 



13j | that all initial data decaying at infinity faster 



than Ae x * x , with A* = 2/VD, will develop travelling waves with velocity c = 1 + 2yD. If 
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n p >0 



FIG. 1: A planar front is perturbed by a displacement of characteristic size 5. The boundary layer 
has a characteristic size y/~D as explained in the text. In this boundary layer, the electron density 
n e goes from 1 to and the positive ion density n p decays to 0. The front moves to the right along 
the positive x axis. 

D<1, the profiles for n p and E will vary very little from the profiles with D = 0. On the 
other hand, n e will develop a boundary layer at the front smoothing the jump from n e = 1 
to n e = 0. Approximating at the boundary layer n p — 0, E — — 1, we obtain the equation 

dn e d 2 n e 

-2— - — = n e (l-n e ) (8) 

where x = ( x ~ (1 + 2v^D)t)/v^D, together with the matching conditions n e {— oo) = 1 and 
n e (+oo) = 0. Eq. (jHJ) is the well known equation for the travelling waves of Fisher's equation. 
It appears in the context of mathematical biology 14] and is known to have solutions subject 
to our matching conditions. This means that we have a boundary layer of width \J~D at 
£ = in which Eq. (jHJ) gives the solution for the electron density n e . Before this layer, we 
have n e m 1, and after the layer, n e w 0. When D = 0, this is the shock front of Eq. (|7|). In 
what respect to n p at the boundary layer, at first order in D one obtains 

/>oo 

n P (x) = -v^D / n e (x)dx, (9) 
J% 

so that drip/dx is 0(1) at the boundary layer. 

Now we make a perturbation in the transversal direction y. We introduce a new system of 
coordinates in the form t = t,y = y,x = x — 5(p(x,y,t) so that, at t = 0, n£\x) , n p ^ (x) , and 
E^ (x) correspond to the profiles of the travelling wave computed in the previous paragraph, 
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and 5 is a sufficiently small parameter (see Fig. Q). By doing this, we follow a strategy 
analogous to the one used in Rubinstein et al. |8[ to deduce the asymptotic approximation 
of Allen-Cahn equation by mean curvature flow. 

We introduce the perturbed electric field and electron and ion densities as 

E = £7<°> u x + 5 u, + E^u y ) , (10) 
n e = n® + 8 2 nW, (11) 
n p = n® + 6n§-\ (12) 

and we will select the function (p(x, y, t) and the 0(5) corrections E x , Ey ,rip in such way 
that Eqs. ([10| ) -([12 p represent correct series expansions. Specifically, the correction to the 
travelling wave profiles for n e will be 0(5 2 ). 

We insert these expressions into Eqs. (0)-©- Then we impose that 0(5°) terms and 
0(5 1 ) terms vanish. The solution of the equations at order 0(5°) is the travelling wave 
found previously, so that rie(x~,t) = rie(x — t) is given by Fisher's equation (JKJ), nf\x — t) 
is given by Eq. Q, and E^°'{x — t) is the solution of Poisson's equation corresponding to 
these particle densities. 

At order 0(5^, after some manipulations, and taking into account that d/dx = 0(D~ l l 2 ) 
and dn^ /dx = 0(1) at the boundary layer, and that D is a small parameter, the equation 
for the perturbed ion density rip decouples from the equations for the perturbed electron 
density and the perturbed electric field (Ex , Ey). The system at 0(5 X ) order is then 
given by the evolution equations 

fUff-g W >- n f)=o. , 14 , 

Observe that the system (|TB |) — ffT^ j) simplifies if one assumes that ip is independent of x. 
This is a valid assumption at least for a short period of time (the one in which stability is 
analyzed) if one assumes ip independent of x at t — 0. 

It is more convenient to formulate Eq. (fTf|) in terms of the electric potential. We note 
that the total electric field has to be irrotational since the magnetic field is negligible. So 
we will assume that E = — VV^, where V is an electric potential that can be written as 
V(x, y) = V^(x) + 5V^(x,y). At order 0(5°), Poisson's equation implies that V^(x) 



is an electric potential associated to the electric field E^(x). At order 0(5 1 ), Poisson's 
equation implies that satisfies 

-V?^)V«(x, y) = -^ m(p —^l _ 2 + n« (15) 

with the condition of decaying at \x\ — > oo. Again, by noting d/dx = 0{D^ 1 / 2 ), we can 
neglect in Eq. f)15|) . and solve the resulting equation by taking Fourier transform in y 
to find the following value for the Fourier transform of Ex , 

^ k)= .m&J 2 ^- m -^° . tie, 

2 { j^e x + j^eW*, for x < 

The front is at a neighbourhood of 0(D 1 ^ 2 ) width around x — 0. While <C D^ 1 ^ 2 , the 
exponentials in Eq. (|16j) can be neglected in this region, and we can write 

^(i )(0 k) = _ (17) 

i.e. a field independent of x. 

Assuming that tp> does not depend on x, we can write Eq. (fT3|) in the form 



The third term at the right hand side of Eq. (fTHjl is the contribution of the electron diffusion 
to the evolution of the front and is also its linearized mean curvature. Taking Fourier 
transform of Eq. (fT%|) in y , and using Eq. (fTTjl , we find 

dip{k) \k\tp(k) 



dt 2(l + |fc|) 
Now, let us write the following ansatz for ip, 



+ D\k\ 2 0{k) =0. (19) 



<p(k,t) = e st (j)(k). (20) 



Introducing this expression into Eq. (JTHJ), we obtain the relation 

lifel 



D\k\ 2 , (21) 



2(1 + |*|) 

that gives the dispersion curve of the transversal perturbations of the planar negative ion- 
ization front explicitly in terms of the parameter D = De/E^ (see Fig. |2J). From this result 
we can obtain some important consequences on the branching of streamers: 




FIG. 2: The dispersion relation given by Eq. 1)21 Jl for D = 0.1 and D = 0.05. For the first curve, 
the position of the maximum (k max ) and the value where the curve changes sign (fep) are shown. 
The physical roles of these values are discussed at the text. 

(i) There exists a maximum of s(\k\) that selects the wavelength of the perturbation. 
When D is a small parameter, this maximum is approximately located at 



1/3 



(22) 



Notice that k max is 0(D -1 / 3 ), so that \k\~x can be safely approximated by zero in the bound- 
ary layer. This justifies the assumption, done previously in this work, that E x l \ and hence 
<p, are independent of x at this order. The value of k max corresponds to a typical spacing 
between fingers given by 



A, 



2tt 



k 1 1 
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£ 



1/3 



(23) 



This is an equation in dimensionless units. If we introduce typical scales for nitrogen (see the 
paragraph below Eq. (j3J), then we can obtain the typical distance between two consecutive 
branches in a negative streamer discharge. This distance will decrease when the gas pressure 
or the initial electric field in the gas increases. 

(ii) There exists an stability threshold, i.e. a value of k delimiting the stability region, 
which corresponds to the nontrivial zero of the function s(k), 



k 



2D P 



The wavelength A associated to this wave number is 

2tt 



An 



ko - 8 - 9 v£ 



(24) 



(25) 



Although the predictions are made for negative planar fronts, they agree with the observed 
fact that, for positive discharges, the number of streamers increases with the electric field 
and the pressure. Those effects are accounted by the expression (j22J). 

We can now provide a qualitative mechanism of streamer branching. When the radius of 
a streamer becomes larger than the critical length given by Eq. (|25|). the streamer becomes 
unstable and branching develops. The electric field which should be taken into Eq. (|23J) in 
case of inhomogeneous electric discharges is the local field at the front of the streamer. In 
order to test the predictions, experimental evidence should be provided in the range where 
the approximations of large electric fields and small diffusion coefficient are valid. 

To conclude, we have derived analytically the characteristic length of branching for planar 
negative ionization fronts. At the same time, this prediction can be considered as a test for 
the validity of the minimal deterministic model on which this calculation is based. 

We thank Robert Deegan for useful discussions on this work. This paper has been par- 
tially supported by the Spanish Ministry of Science and Technology grant BFM2002-02042, 
and by the Universidad Rey Juan Carlos grant PPR-2004-38. 



[1] N. Liu and V. P. Pasko, J. Geophys. Res. 109, A04301 (2004). 
[2] C. Montjin, private communication. 

[3] M. Arrayas, U. Ebert and W. Hundsdorfer, Phys. Rev. Lett. 88, 174502 (2002). 
[4] V. P. Pasko, M. A. Stanley, J. D. Mathews, U. S. Inan and T. G. Wood, Nature 416, 152 
(2002). 

[5] M. Arrayas and U. Ebert, Phys. Rev. E 69, 036214 (2004). 

[6] A. N. Kolmogorov, I. G. Petrovskii and N. S. Piskunov, Study of the diffusion equation with 
growth of the quantity of matter and its application to a biology problem, in Selected Works 
of A. N. Kolmogorov (Kluwer Academic Publishers, Amsterdam, 1991). 

[7] S. M. Allen and J. W. Cahn, Acta. Metal. 27, 1085 (1979). 

[8] J. Rubinstein, P. Sternberg and J. B. Keller, SIAM J. Appl. Math. 49, 116 (1989). 

[9] J. W. Cahn and J. E. Hilliard, J. Chem. Phys. 28, 258 (1958). 
[10] N. D. Alikakos, P. W. Bates and X. Chen, Arch. Rat. Mech. Anal. 128, 165 (1994). 
[11] E. N. Dancer, D. Hilhorst, M. Mimura, and L. A. Peletier, Eur. J. Appl. Math. 10, 97 (1999). 



9 



[12] M. Arrayas, M. A. Fontelos and J. L. Trueba, http://xxx.lanl.gov/abs/physics/0504005, 
[13] U. Ebert, W. van Saarloos and C. Caroli, Phys. Rev. Lett. 77, 4178 (1996); and Phys. Rev. 
E 55, 1530 (1997). 

[14] J. D. Murray, Mathematical Biology (Springer- Verlag, New York, 1990). 



10 



